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INTRODUCTION 


The  present  report  originated  from  a  specific  question: 
can  the  effectiveness  of  a  jet  propulsion  motor  be  increased 
by  controlling  the  mode  of  detonation  or  combustion?  Methods 
previously  developed  for  gas  dynamical  problems  by  the  Applied 
Mathematics  Group  of  New  York  University  provide  a  basis  for 
making  a  theoretical  contribution  by  computing  the  pressure 
as  a  function  of  the  time  under  highly  idealized  assumptions, 
thereby  comparing  the  effectiveness  of  different  modes  of  reac- 
tion.  Such  computations  are  presented  in  the  present  memorandum. 

However,  the  value  of  the  study  should  not  be  judged 
solely  in  relation  to  the  original  specific  problem.   Methods 
and  general  results  may  well  be  useful  for  many  other  problems 
of  detonation  and  combustion. 

The  work  was  carried  out  by  Dr.  Charles  De  Prima  and 
Miss  Anneli  Leopold  under  the  guidance  of  Dr.  K.  0.  Priedrichs. 
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On  the  Effectiveness  of  Various  Modes  of 


Detonation  or  Combustion 


1.  Problem*   In  connection  with  the  method  of  cascading  In 
Jetmotors  proposed  by  Dr.  P.  Zwlcky  from  the  Aerojet  Co. 
[1],[2]  we  have  Investigated  the  effectiveness  of  various 
modes  of  detonation  or  combustion  on  the  basis  of  a  simpli- 
fied model. 

We  have  assumed  that  a  straight  tube,  extending  to 
Infinity  In  one  direction  and  closed  at  the  other  end,  contains 
explosive  gas  In  a  section  adjacent  to  the  closed  end  and  air 
In  the  remaining  part.   Then  we  have  considered  the  following 
four  modes  of  detonation  and  two  modes  of  combustion. 

DI:   The  explosive  gas  Is  Ignited  at  the  closed  end  in  such  a 
way  that  a  detonation  wave  is  produced  which  sweeps  over  the 
section  of  explosive  gas.   (See  Figure  1). 


Figure  1 


DII:   The  explosive  gas  is  ignited  at  the  interface  between 
air  and  explosive  in  such  a  way  that  a  detonation  wave  is 
produced  which  moves  toward  the  closed  end.  (See  Figure  2). 


Figure  2 


Dili:  The  whole  section  of  explosive  is  ignited  instantaneous- 
ly so  that  a  constant  volume  detonation  results. 

DIV:  The  explosive  gas  is  first  adiabatically  compressed  by 
the  action  of  a  piston,  then  ignited  Instantaneously.  Here- 
after the  exploded  gas  is  allowed  to  expand  adiabatically  but 
is  forced  to  work  against  the  piston  until  it  has  delivered  the 
work  it  has  received  during  the  pre-compression  phase;  only  then 
the  piston  is  removed  and  the  exploded  gas  may  expand  pushing 
out  the  air  ahead  of  it.   (See  Figure  3). 

CI:  Ignition  at  the  closed  end  produces  an  ordinary  combustion 
wave  moving  away  from  the  closed  end.   (Pig.l) 

CII:  Ignition  at  the  interface  produces  an  ordinary  combustion 
wave  moving  toward  the  closed  end.   (Fig.  2) 

To  measure  the  effectiveness   of  these  modes  of  opera- 
tion we  have  calculated  the  impulse  per  unit  area  delivered  to 
the  closed  end  i.e.  the  Integral  of  the  excess  pressure  with 


We  use  the  somewhat  vague  term  effectiveness  instead  of 
efficiency  to  avoid  confusion  since  the  latter  terra  is  general- 
ly used  for  the  ratio  of  useful  energy  output  to  energy  input. 


respect  to  the  time, 

I  =  I(t)  =/  [p-pjdt, 
o 

p  being  the  pressure  exerted  by  the  gas  against  the  closed  end 
and  Pq  being  the  atmospheric  pressure. 

Our  calculations  extend  only  over  a  short  Interval  of 
time;  but  this  Interval  Is  large  enough  to  cover  the  decisive 
part  of  the  process:  the  sudden  rise  and  the  quick  drop  of  the 
pressure.   After  the  time  up  to  which  our  calculations  extend, 
the  pressure  decreases  very  slowly.   In  all  modes  of  operation 
considered  we  have  assumed  that  Initially  the  same  amount  of 
combustible  gas  is  available.   Tghen  we  say  that  one  mode  of 
operation  Is  more  effective  than  another  we  mean  merely  that 
this  mode  produces  a  higher  Impulse  than  the  other  during  the 
time  interval  considered. 

The  reason  for  considering  impulse  rather  than  energy  as 
a  measure  of  effectiveness  is  this:   For  a  Jet  motor  which 
operates  in  a  periodic  process,  the  average  thrust  over  a  period 
may  be  considered  a  measure  for  effectiveness.   This  average 
thrust  is  given  by  the  expression 

^''(P-Po)^*^ 
■■■  o 

if  T  is  the  period.   The  main  contribution  to  the  value  of  the 
integral  here  is  due  to  a  short  interval  of  time  after  the  ex- 
plosion, during  which  the  pressure  rises  suddenly  and  then  drops 
again.   Thus,  when  one  compares  two  processes  operating  with  the 
same  period  T,  one  may  assume  that  the  more  effective  process  is 
that  one  having  the  higher  value  of  the  impulse 


I  =/  (P-Po)dt 


o 
over  the  first  decisive  interval  of  time. 
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Naturally,  for  the  discussion  of  the  operation  of  a 
Jet-propelled  device,  completely  different  notions  of  effective- 
ness". In  particular  an  over-all  l^t/tgfi^^7  b©  needed.   Such  an 
over-all  efficiency  depends  on  various  factors  not  considered 
here.   Our  sole  Intention  In  the  present  report  Is  to  furnish 
only  some  Information  about  one  of  the  factors  to  be  taken  Into 
account  when  such  an  over-all  efficiency  Is  Investigated. 

In  our  calculations  we  hay.e  made  particular  assumptions 
about  the  nature  of  the  explosive  gas  which  we  shall  explain  In 
detail  later  on.   The  Initial  state  was  assumed  at  atmospheric 
pressure  and  room  temperature. 

2.  Results.   The  results  In  general  terms  are  the  following: 

The  Impulse  resulting  from  mode  DII  exceeds,  after  a 
short  time,  that  resulting  from  mode  DI  and  remains  greater  by 
approximately  a  constant  value  for  some  time.   In  this  sense 
It  Is  more  effective  to  let  the  detonation  proceed  toward  the 
wall  than  away  from  It.   However,  the  difference  Is  not  very 
great. 

To  express  the  result  numerically  In  dlmenslonless  form 
we  must  Introduce  a  "unit  Impulse"  I  ;  to  this  end  we  Introduce 
the  atmospheric  pressure  p  as  unit  pressure,  the  length  1   of 
the  explosive  section  as  unit  length,  and  as  unit  time  the  time 
tp  =  IqAo  which  a  sound  wave  in  atmospheric  air  would  take  to 
travel  across  the  unit  length.   If,  for  example,  the  unit  length 
1q  =  34.7  cm.  =  1.14  ft.,  then  the  unit  time  is  just  one  milli- 
second.  The  unit  Impulse  I   is  then  the  Impulse  p.u.  area  which 
atmospheric  pressure  would  produce  in  unit  time, 

^o  =  Po^o  =  Po^oAo  • 

We  find,  for  the  particular  case  considered,  described 
later  on,  that  at  the  time  t  =  tQ  the  impulse  from  mode  DI  has 

*  For  various  notions  of  efficiency  in  connection  with  jet- 
propelled  devices  sec  Zwlcky  [2], [3]. 


attained  the  value  I  =  2.91  ,  and  that  from  mode  DII  the 
value  I  =  5.21^;  from  there  on  the  values  of  I  increase  only 
slowly  and  the  difference  of  the  two  values  of  the  impulse 
remains  approximately  constant,  at  about  the  value  21   for 
some  time.   More  details  can  be  read  off  from  Table  1. 

The  speed  of  the  detonation  waves  in  the  modes  Dl  and 
DII  is  the  "natural"  one  determined  on  the  basis  of  the  Chapman- 
Jouguet  hypothesis  (see  Manual  [4]  pp.  129-133  and  Supplement 
[5],  Lewis  and  v.  Elbe  [6]  p. 235  and  Jost  [7]  p.  157).   Other 
modes  of  detonation  would  result  from  an  artificially  propagated 
detonation.   The  simplest  such  artificial  mode  is  the  constant 
volume  detonation,  Dili,  initiated  by  instantaneous  ignition. 
The  value  of  the  impulse  in  this  case  remains  between  those  of 
the  modes  Dl  and  DII  for  some  time.   We  could  not  ascertain 
whether  or  not  it  eventually  gets  ahead  of  both  and  increases 
at  a  faster  rate  than  those  of  the  other  two  modes.   An  investi- 
gation of  mode  Dili  is  justified  mainly  because  the  calculations 
for  its  determination  are  relatively  simpler;  not  because  in- 
stantaneous detonation  might  be  produced  in  actual  performance. 
We  were  Interested  in  particular  In  finding  out  to  which  degree 
of  accuracy  the  results  of  any  mode  of  detonation  could  be  re- 
placed by  those  of  a  constant  volume  detonation;  and  we  have 
found  that  this  can  be  done  with  good  accuracy  (see  Tables  1  and 
3,  representing  Impulse  and  pressure  as  functions  of  time  for  the 
three  modes  Dl,  DII,  Dili.) 

Still  a  different  mode  of  detonation  had  been  suggested: 
First  a  section  of  the  explosive  adjacent  to  the  interface  la 
detonated  Instantaneously;  a  second  section  is  detonated  after  it 
has  been  compressed  by  the  shock  wave  sent  out  from  the  first 
detonation.   Successively,  one  section  after  another  is  detonated 
after  it  has  undergone  one  shock  compression  after  another  result- 
ing from  the  previous  detonations.   Due  to  the  high  pre -compress ion 
of  the  last  sections,  a  high  increase  of  efficiency  could  be 


expected.  It  seems  most  likely,  however,  that  the  Impulse 
resulting  from  such  a  process  lies  between  those  resulting 
from  modes  DII  and  Dili. 

To  obtain  the  benefit  from  pre -compress ion  it  seems 
necessary  to  produce  the  pre -compress ion  by  external  means. 
Therefore,  the  idealized  mode  DIV  has  been  investigated. ' 
Various  models  may  be  imagined  to  produce  it.   In  any  case, 
there  are  the  following  four  preliminary  steps  in  this  ideal- 
ized mode,  as  indicated  in  Figure  3.   Initially,  the  explosive 
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has  atmospheric  pressure  Pq,  and  volume  V^.   A  piston  placed 
at  the  Interface  between  explosive  and  air  then  compresses  the 
explosive  adiabatically  to  a  prescribed  vol\ime  v'  and  to  a 
pressure  p'  (see  Pig.  3,  A  and  B).   The  piston  is  to  come  to  rest 
when  this  volume  is  reached.   At  this  instant  the  explosive  id 
ignited  instantaneously,  thus  acquiring  the  pressure  p|  but  still 
occupying  the  volume  v|  =  v'  (see  Pig.  3C).   The  piston  then 
moves  backwards  and  the  exploded  gas  expands,  adiabatically  by 
assumption  (see  Pig.  3D).   As  soon  as  the  exploded  gas  has,  in 
pushing  back  the  piston,  done  as  much  work  as  the  piston  had 
done  In  compressing  the  gas  from  p   to  p'   the  piston  is  removed. 


*  It  was  suggested  to  us  by  Mr.  E.  S.  Roberts. 


It  is  easily  seen  that  the  total  Impulse  Imparted  to  the  tube 
during  this  preliminary  period  of  the  process  Is  zero.   Prom 
the  time  on  at  which  the  piston  was  removed  (we  shall  refer 
to  this  time  as  the  "effective  Initial  time"),  the  exploded 
gas  interacts  with  the  air  in  the  same  way  as  after  a  constant 
volume  detonation  without  pre-compression.   We  have  measured 
the  impulse  from  this  effective  initial  time  on.   Pressure  and 
volume  at  this  time  are  denoted  by  p,  and  V,  respectively. 

This  Impulse  was  calculated  for  various  degrees  of  pre- 
compression  viz.  for  the  ratios 

V^^/V^  =  1.00,  6.03,  103.2 

of   original  volume   V     to  pre-compreased   volume   V      (see  Fig.    3, 
AjB).      The   resultant   Impulses   are   shown   in  Table   2,    the  resultant 
pressures    in  Table  4.      It   is    seen  that  with  greater  pre-compres- 
sion  Indeed  a  greater   Impulse  results;   but   the   Increase   is   not 
very  large.      While   for  case   V  /v'    =  1,    the   impulse  attains    the 
value    I   =  4.75  I^   at    t   =  t^   and   I    =  7.6   I^   at    t    =  3t^,    for   V/v'  = 

O  O  0  O'        O'   O 

6.03,  the  impulse  is  I  =  7.55  I   at  t  =  t  and  I  =  9.7  I   at 

0  0  0 

t  =  3  t  •   The  difference  thus  changed  little  during  this  time 
and  remains  approximately  constant  at  the  value  2.1  I   from  then 
on.   It  is  perhaps  surprising  that  this  difference  is  relatively 
low,  and  that  the  impulse  increases  only  little  if  the  pre-compres- 
sion is  further  increased,  even  to  values  of  the  compression  ratio 
as  high  as  V  /v'  =  103.2  as  seen  from  Table  4.   In  this  comparison 
the  duration  of  the  preliminary  (from  the  beginning  of  the  pre- 
compression  to  the  effective  initial  time)  phase  has  been  ignored. 
•This  duration  depends  on  the  mechanism  employed,  but  one  may 
assume  that  it  will  Increase  with  the  degree  of  compression  V/V  . 
Therefore,  if  this  duration  is  long,  the  Impulse  will  become 
larger  than  that  from  a  constant  volume  detonation  only  after  an 
accordingly  long  time. 
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It  was  natural  to  expect  that  pre-ccmpresslon  would  in- 
crease the  effectiveness  considerably  since  this  is  the  case  in 
a  combustion  motor.   There  is,  however,  a  decisive  difference 
between  the  process  in  a  combustion  motor  and  the  process  con- 
sidered here.   In  a  combustion  motor,  the  exploded  gas  is  allowed 
to  move  the  piston  back  to  the  original  position,  thus  perform- 
ing excess  work  from  the  "effective  initial  time"  on.   This  ex- 
cess work,  which  is  considerable  for  a  high  pre-ccanpression,  is 
used  to  drive  an  augmentation  device  which  moves  large  masses; 
such  as  air  through  an  airplane  propeller,  water  through  a  ship 
screw,  or  the  "whole  earth"  by  a  land  vehicle",  and  it  is  only 
due  to  this  augmentation  that  the  thrust  is  increased.   In  the 
present  process,  however,  in  which  the  piston  is  removed  at  the 
effective  initial  time,  the  energy  available  for  the  excess  work 
is  only  used  to  accellerate  and  heat  the  exhaust  gases  so  that 
no  great  benefit  is  derived  from  the  pre-compression.   Thus  it 
seems  that  one  should  not  expect  that  pre-compression  in  a  jet 
motor  could  step  up  the  impulse  to  the  same  degree  as  in  a  com- 
bustion motor  unless  considerable  augmentation  is  used. 

This  somewhat  negative  statement  should,  however,  not  be 
over-rated.   Our  calculation  concerns  only  the  first  phase  of 
the  operation.   The  situation  is  perhaps  different  if  the  whole 
cycle  is  considered.   The  following  argument  might  be  offered 
in  support  of  this  idea.   A  rocket  motor  delivering  (at  least 
over  a  considerable  length  of  time),  a  steady  exhaust  flow 
through  a  nozzle  with  atmospheric  pressure  at  the  exhaust,  pro- 
duces a  tbjrust  given  by  P  =  mv,  m  being  the  mass  flux  per  unit 
time  and  v  the  exhaust  gas  velocity;  the  exhaust  flux  of  kinetic 
energy  is  given  by  —  mv  .   Therefore  any  device  which  increases 
the  flux  of  kinetic  energy  automatically  increases  the  thrust 
and  hence  the  impulse. For  a  motor  with  periodic  operation,  the  im- 

pulse  over  a  cycle  of  period  T  is  /  mvdt  and  the  kinetic  energy 

o 


«•  Cf.  P.  Zwicky,  [2],  page  31, 


1   T.  2 
ejected  during  a  cycle  Is  -^  /  mv  dt.   It  la  no  longer  evident 

'^  o 
that  Increase  of  the  ejected  kinetic  energy  Implies  Increase 

of  Impulse,  but  It  Is  still  likely  to  a  certain  degree.   And  since 

a  pre -compress ion  device  outlined  above  Increases  probably  the 

kinetic  energy  ejected  during  a  cycle,  it  probably  also  increases 

the  Impulse. 

A  completely  different  process  results  if  the  ignition  does 
not  produce  a  detonation  but  only  an  ordinary  combustion  process. 

We  have  made  the  assumption  that  the  width  of  the  zone 
across  which  combustion  takes  place  is  negligibly  small,  in  other 
words,  that  combustion  takes  place  across  a  discontinuity  surface; 
such  combustion  processes  will  also  be  called  "deflagrations."  In 
most  actual  cases,  the  combustion  zone  does  not  only  have  a  notica- 
able  width,  but  in  general  differs  considerably  from  a  plane.   Fre- 
quently, the  flow  involving  combustion  is  turbulent  and  not  steady. 
Nevertheless,  our  assumption  can  be  employed  to  describe  approxi- 
mately the  situation  even  in  such  cases  (although  the  approximation 
is  certainly  crude)  provided  that  burning  takes  place  in  a  finite 
zone  whose  width  remains  nearly  constant  in  time  and  at  whose  two 
sides  the  quantities  are  almost  constant  over  the  cross -section. 
For,  in  that  case  the  conservation  laws  give  the  same  relation  for 
the  quantities  at  both  sides  of  the  zone  as  for  both  sides  of  a 
discontinuity  surface.   A  sufficiently  large  neighborhood  of  a 
curved  burning  surface  or  a  region  of  turbulent  burning  may,  to  a 
certain  degree,  be  considered  such  a  zone.   (Of  course,  this  argu- 
ment assumes  that  losses  due  to  viscous  friction  and  due  to  heat 
conduction  through  the  wall  can  be  ioinored.  ) 

The  decisive  difference  between  the  propagation  of  a  plane 
burning  front,  of  a  curved  burning  front,  and  of  a  laminar  or 
turbulent  narrow  burning  zone  whose  width  remains  constant  in  time 
lies  in  the  speed  of  propagation.  The  (normal)  burning  speed, c«, 
of  a  flame,  i.e.  the  speed  of  the  front  relative  to  the  unburnt 
gas  ahead  of  it  dejiends  predominantly  on  the  heat  conductivity  ot 
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the  gases.   Its  determination  involves  principles  quite  differ- 
ent from  those  determining  the  natural  detonation  speed;  (cf.  Lewis, 
V.  Elbe  [6],Jost  [7]p.63ff  and  a  forthcoming  report  of  the  NYTJ-AMG). 
For  a  plane  burning  front,  the  value  of  the  burning  speed  is  very 
low,  of  the  order  of  magnitude  of  .5  to  3  m/sec  For  curved  burn- 
ing fronts,  the  local  burning  speed  is  the  same  as  for  a  plane 
front;  but  when  we  describe  the  propagation  of  a  curved  flame 
front  approximately  by  that  of  a  plane  burning  zone,  the  "effective" 
burning  speed,  c-,  corresponding  to  such  a  description  is  greater. 
The  same  is  true,  to  a  higher  degree,  for  turbulent  burning.  (Cf» 
Jost  [7]p.  63ff,  89ff).   (Jost  mentions  in  his  book  on  p. 61  that 
flame  speeds  may  be  of  the  order  of  magnitude  of  several  hundred 
m/sec).   In  our  calculations  we  have  assumed  various  values  for 
the  effective  burning  speed,  viz. 


c^/c^  =  .01,  =  .05,  =  .15 


or 


„   =   3.47  m/sec,    =   17.35  m/sec,    =   52.05  m/sec. 


c-  = 


As  indicated  in  the  introduction  we  have  considered  two 
modes  of  burning:  CI,  in  which  the  flame  begins  at  the  closed  end 
and  CII,  in  which  the  flame  begins  at  the  interface  and  moves  to- 
ward the  closed  end.   It  is  clear  that  in  the  case  CII  the  pressure 
at  the  closed  end  starts  to  build  up  only  after  the  pre -compress ion 
wave  sent  out  by  the  flame  has  hit  the  closed  end;  therefore  the 
impulse  resulting  from  mode  CI  will  for  some  time  remain  greater 
than  that  resulting  from  mode  CII;  eventually,  however,  the  im- 
pulse resulting  /rom  mode  CII  will  be  greater  than  that  from  mode 

CI.   This  is  seen  on  Table  1  for  c^/c   =  .15.   If  the  effective 

f '  o 

burning  velocity  c^  is  low,  c^/c^  =  .01  or  =  .05,  the  process  is  so 
slow  that  no  appreciable  impulse  will  have  built  up  during  the 
time  t  =  3t^  considered  by  us.   It  is  most  remarkable  that  for  the 
higher  burning  speed,  c^  =  -ISc^,  the  impulse  from  mode  CII  assumes 
the  value  I  =  5.251^  at  t  =  ^.-^t^  and,  hence,  is  larger  than  that 
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from  the  detonation  mode  DI,  at  that  time.   Its  value  will  most 
likely  never  exceed  that  resulting  froih  mode  DII;  It  Is  at  least 
certain  that  the  rate  of  increase  of  the  Impulse,  I.e.  the  ex- 
cess pressure,  will  drop  slowly  from  t  =  1.23t   on  and  consider- 
ably from  about  t  =  3t   on. 

o 

It  is  of  interest  to  note  that  the  flame  speeds  U-  in  these 
processes,  i.e.  the  speeds  of  the  flame  relative  to  the  tube  are 
much  larger  than  the  effective  burning  speeds.   The  values,  which 
are  different  for  the  two  modes  CI  and  CII  are  indicated  in  the 
Table  below. 

c^/Cq  =   .01        .05      .15 
or  in  n/sec,        c^  =  3.47    ,   17.35  ,   52.05  . 

CI  ^f/^o  ^    *^^  '^^  '^'^ 

or,  in  n/sec,       U^  =  27.7    ,   125.    ,  336. 

U^/c^  =    .04        .22      .70 


CII 


or,  in  ny'sec,       U~  =  15.2    ,   76.5   ,  250. 


Thus  in  the  case  c-  =  'ISc   the  flame  moves  with  approximately 
sound  speed  in  mode  CI. 

3.  Properties  of  the  Gases.   It  was  not  our  intention  to  assume 
a  specific  well-defined  explosive  and  make  precise  predictions 
which  should  be  verified  experimentally.   Our  only  intention  was 
to  obtain  typical  results.   We  have  assumed  the  explosive  to  be  an 
ideal  gas  whose  temperature  T  is  given  by 

T  =  m"-^r  e 

where  M  is  the  molecular  weight,  R  the  gas  constant  and  9  the 
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"reduced  temperature" 

9  =  pr  =  pp"^   , 

p  being  the  pressure,  r  the  specific  volume  (volume  per  unit  mass), 
P  the  density  (mass  per  unit  volume).   We  have  further  assumed 
that  the  internal  energy  per  unit  mass,  e,  depends  only  on  the  re- 
duced temperature  9,  and  not  on  the  pressure;  assuming  the  molecu- 
lar weight  to  be  constant,  or  to  depend  only  on  9  ,  this  means 
that  the  internal  energy  depends  only  on  the  temperature.   Denoting 
the  energy  of  formation  per  unit  mass  (at  absolute  zero  temperature) 
by  g  we  have  for  the  total  energy  per  unit  mass 

E  =  e  +  g  . 

For  air  we  simply  set  g  =  0  and 

^  =rV*  ,  «^  =  1.4 . 

For   the    initial  temperature   T     =  500°K  we  have   9     =  .0206   kcal/gm. 
For   the  unburnt   gas   we  have  assumed   (setting  g      =0) 

with     ^      =   1.43.      (We  have   chosen    <f     =   1.43    ir^tead   of   1.4  because 
the  unburnt   gas    is   a  mixture   of  air  with  a   small   fraction  of   oil 
droplets   which  affect   the   weight  somewhat  but   not    the   volume). 
For  the  burnt  gas   we  have   adopted  the   relation 

E   =   -   .679   +   2.59   +   4.569^   +    19.869*^    -   36.829'^ 

where  9  and  E  refer  to  kcal/gm  as  units.   The  energy  of  formation 
of  the  burnt  gas  (at  the  absolute  temperature  zero)  is 

g^  =  -  .679  . 
Therefore,  the  "liberated  energy"  is 
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The   ratio  of  specific  heats    is    to  be   found   from 

_1_  =  H  =  2.5   +   9.129   +   59.580^    -   147.289^    . 
V-1        ^ 

For  T  =  300°,  900°,  and  3000°  K  we  have   ^=  1.4,  1.51,  and 
1.18,  respectively,  see  Table  6  where  E  and  ^  are  plotted 
as  functions  of  0  or  T.   In  the  relationship  between  the 
temperature  T  and  the  reduced  temperature  0  we  have  assumed 
M  =  29  throughout  ignoring  the  slight  differences  of  the 
values  of  M  for  air,  unburnt  gas,  and  burnt  gas  at  different 
temperatures. 

The  particular  expression  for  E  as  function  of  0  was 
obtained  by  assuming  the  explosive  to  consist  of  a  mixture  of 
air  and  hydrocarbon  such  that  the  numbers  of  carbon  atoms,  hy- 
drogen atoms,  oxygen  molecules,  and  nitrogen  molecules  were 
in  the  ratio  1:1.78:1.45:5.50.   In  determining  the  energy 
expression  the  change  in  composition  Involving  some  dissocia- 
tions was  taken  into  account;   the  deviation  of  the  expression 
for  E  from  linearity  in  0  is  mainly  due  to  the  changes  of  the 
specific  heats.   We  should  like  to  emphasize  once  more,  that 
our  results  do  not  pretend  to  represent  the  behavior  of  a 
particular  substance  accurately;   they  are  only  intended  to  be 
typical.   As  a  matter  of  fact,  calculations  have  been  performed 
previously  under  the  assumption  that  air,  unburnt,  and  burnt 
gas  have  the  same  constant  specific  heats  with  the  ratio 
0     =   1.334,  the  liberated  energy  -  g,  being  approximately 
•75  kcal/gm.   The  results  were  essentially  the  same  as  in  the 
case  here  reported;  except  that  the  values  of  the  Impulse  are 
somewhat  higher  and  that  the  value  of  the  impulse  for  case  Dili 
is  during  some  time  nearer  to  that  for  case  DI  rather  than  to 
that  for  case  DII. 

4.  Detailed  Description  of  the  Processes.  We  proceed  to  describe 
the  motion  of  the  p;ases  and  the  variation  of  pressure  and  teyiper- 
ature  resulting  from  the  various  modes  of  detonation  and  combus- 
tion.  The  flow  processes  will  be  described  in  terms  of  "waves". 
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Detonation,  combustion,  and  shock  waves  are  nothing  but 
moving  detonation,  combustion,  and  shock  fronts  Implying 
a  discontinuity.   Rarefaction  and  compression  waves  are 
moving  zones  across  which  pressure,  temperature,  and 
velocity  vary.   That  such  a  description  is  suitable  for 
actually  occuj'rlng  processes  ma^/  be  seen  from  Schlieren- 
photographs  made  by  Payman  and  his  collaborators. (See  [8]). 
We  first  consider  the  process  resulting  from  a  constant  volume 
detonation,  mode  Dili.   After  the  detonation  the  exploded  gas 
has  Instantaneously  acquired  a  pressure  p  =  9.56  atm.  and  a 
temperature  T  =  2800°K.   Since  the  pressure  difference  between 
9.56  atm.  and  1  atm.  at  the  interface  between  exploded  gas  and 
air  cannot  be  maintained,  waves  are  sent  out  from  this  inter- 
face which  produce  an  adjustment.   A  shock  wave  travels  into 
the  quiet  air  raising  its  pressure  to  5.91  atm.,  its  tempera- 
ture to  542^  and  accelerating  it  to  a  velocity  u  =  1.4c  ; 
here  c   =  347m/sec  =  1.14  ft/sec  Is  the  sound  speed  in  air  of 
atmospheric  pressure  and  300  K  temperature.   A  rarefaction 
wave  travels  backwards  into  the  exploded  gas,  reducing  the 
pressure  to  5.91  atm.,  the  temperature  to  2550  a  and  accelerat- 
ing the  exploded  gises  to  a  velocity  u  =  1.4c  .   The  rarefaction 
wave,  in  contrast  to  a  shock  wave,  has  a  finite  width  which  in- 
creases proportional  to  the  time  it  has  traveled.   This  rare- 
faction wave  will  be  reflected  at  the  closed  end;   due  to  this 

reflection,  which  takes  place  in  the  time  interval  between  t=.35t 
*  o 

and  t  =  .67t  ,  the  pressure  at  the  closed  end  is  reduced  from 

9.36  atm.  to  3.25  atm.  and  the  temperature  from  2800°K  to  2300 °K. 

The  reflected  rarefaction  wave  meets  the  interface  between  burnt 

gas  and  air  at  the  time  t  =  .94t  .   At  this  time,  the  interface 

has  moved  out  to  the  distance  2.31  1   from  the  closed  end.   The 

o 

reflected  rarefaction  wave  will  be  partly  transmitted  and  partly 
reflected  at  the  interface.   The  effect  of  this  Interaction  is 
a  considerable  deceleration  of  the  Interface.   In  the  time  Interval 
t  =  •94t   to  1.83t  ,  the  velocity  decreases  from  1.4c   to  .57c  . 
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At  the  same  time  the  temperature  Is  reduced  to  2170  K  In  the 
burnt  gas  and  to  423  K  In  the  air.   The  rarefaction  wave  trans- 
mitted through  the  interface  will  catch  up  with  the  shock  wave 
(beginning  at  t  =  2.1  t   at  a  distance  5.5  1   from  the  wall). 
The  Interaction  between  rarefaction  wave  and  shock  wave  con- 
sists of  a  gradual  diminishing  and  slowing  down  of  the  shock 
wave.   This  process  Is  extremely  slow  and,  probably,  goes  on 
Indefinitely;   the  reflected  wave  resulting  from  this  Inter- 
action Is  negligible.   However,  the  rarefaction  wave  reflected 
from  the  interface  will  again  be  reflected  at  the  closed  end; 
the  resulting  reflected  wave  will  again  be  partly  reflected, 
partly  transmitted  at  the  Interface,  and  so  on.   Three  such 
reflections  at  the  closed  end  were  determined.   The  pressure 
and  the  temperature  at  the  closed  end  will  gradually  decrease 
In  the  course  of  these  reflections;   the  pressure  there  will 
probably  approach  atmospheric  pressure  asymptotically.   This 
decrease  of  pressure  is  shown  In  Table  3. 

The  results  here  described  are  exhibited  In  the  diagram. 
Table  7,  in  which  the  distance  x  from  the  closed  end  is  plotted 

as  absicca  and  the  time  t  as  ordinate.   The  units  are  1  and 

o 

t^  =  1q/Cq  respectively.   The  shock  is  indicated  by  a  heavy  line, 
the  interface  by  a  dotted  line.   To  plot  rarefaction  waves,  we 
have  made  use  of  the  fact  that  the  end  points  of  the  wave  zone 
move  with  sound  velocity  (relative  to  the  gas),  which  of  course 
is  different  at  both  end  points.   The  region  in  the  (x,t) -plane 
covered  by  a  rarefaction  wave  is  therefore  bordered  by  lines 
which  indicate  the  path  a  sound  wave  (i.e.  an  extremely  weak 
shock  or  rarefaction  wave)  would  describe.   Such  lines  we  shall 
refer  to  as  "sound  lines."   It  is  then  natural  to  describe  a 
rarefaction  wave  by  showing  several  sound  lines  in  the  region 
covered  in  addition  to  those  bordering  this  region. 

The  flow  just  described  was  determined  on  the  basis  of 
the  assumption  that  heat  conduction  and  viscosity  can  be  Ignored 
(despite  the  fact  that  both  play  a  role  in  the  creation  of  a 
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shock).      It  appears   to  be  rather   severe   to   ignore  heat   conduc- 
tion in  the  neighborhood   of   the    Interface  between  burnt  gas 
and  air   In  view  of   the   fact   that   the   difference   of  the    tempera- 
tures   on  both  sides    Is   2000 °K  at    first   and  still   over   1600 °K 
after  3   time  units. 

Next  we  describe   the  processes   resulting   from  the  mode  DIV, 
Involving  pre-compresslon.    For  the   case   of   the   pre-compresslon 
ratio  V  /V'    =  6.03,    the   gas   has   acquired  the   pressure   12.06   atm. 
and   the    temperature  600°K;    after   the    explosion  the   pressure    Is 
58.4  atm.    and  the   temperature   2900   .      In  the   subsequent   pre- 
expanslon  leading   to  ".state    (1)"    the  volume   Is    only  Increased 
from  V»^  =  .166   V^^  to  V^  =   .198  V^   (see   Figure   3);    the  pressure 
drops    to  47.32  atm.    and   the    temperature    to  2800 °K.      Incidentally, 
this   latter   temperature   is    independent   of  the  degree   of  pre- 
ccmpression;      the   reason  is    that   the  value   of  the   total  energy 
E  =  e  +  g   in  state    (1)    is   equal   to   that    in  the   original  state    (0) 
since   the   energy   increase   during  pre-compresslon  is   balanced  by 
the    energy   loss   during  pre-expanslon  and  the  total   energy  does 
not   change   in  a   constant  volume   detonation.      When  the   state    (1) 
Is   reached,    at   the   "effective   initial   time,"   the  piston  is   re- 
moved and  a   flow  process    ensues  which   is   similar   to  that  result- 
ing  from  a   constant   volume   detonation  without  pre-compresslon. 
The   only  difference    is    that   the    initial  volume    is   smaller, 
V,    =  .198  V   ,    the   initial  pressure  p,    =  47.37  atmo    is  higher, 
while   the    initial   temperature   is   the   same.      In  consequence    of 
these  differences    the  various   reflections   take  place   faster. 
We  have   determined   three   reflections   at   the   closed   end,    the   last 
being   completed  at    t  =  5.54   t    j      the  pressure  has    dropped  to 
1.4   atm.    at    that    time    (see   Table   8). 


<•  Calculations  made    on   the   basis    of    the    standard  heat    flow 
theory   indicate,    however,    that  after   3   time   units    the  width 
of  the  region  in  which   80  percent    of  the   temperature   difference 
is   attained   is    only   of  the   order   of  magnitude    of   one   cm. 
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In  the  extreme  case  of  the  compression  ratio  V  /V  =103.2, 
the  volume  Is  v'  =  .0097  V  ,  the  pressure  Is  515.8  atm. ,  the 
temperature  1500  K.   As  a  result  of  the  explosion,  the  pressure 
rises  to  1176  atm.  and  the  temperature  to  3420  K.   Through  pre- 
©xpanslon,  the  volume  Increases  to  V,  =  .0425  V  =  (1/23. 5)V^, 
the  pressure  drops  to  200.15  atm.  while  the  temperature  again 
Is  2800°K.   The  reflections  take  place  very  quickly.   The  first 
reflection  at  the  wall  Is  completed  at  t  =  .22  t^.   It  decreases 
the  pressure  to  2.8  atm.,  the  temperature  to  1250  K.   These 
values  remain  there  until  the  time  t  =  .44  t  when  the  next  re- 
flection begins.   We  have  stopped  our  calculations  there  (see 
Table  9).   The  Impulse  as  function  of  the  time  t  up  to  t  =  .44  t^ 
plotted  In  Table  1,  already  shows  clearly  that  the  gain  In  Im- 
pulse Is  slight  as  compared  to  the  tremendous  Increase  In  the  pre- 
compresslon  rate.  (Rough  calculations  have  Indicated  that  the 
rate  of  increase  In  Impulse  as  function  of  the  pre -compress ion 
rate  begins  to  drop  when  its  value  is  near  V  /V  =  6. ) 

The  flow  processes  resulting  when  a  detonation  wave  travels 
either  away  from  the  wall  or  toward  the  wall  are  slightly  more 
complicated.   In  mode  DI,  see  Table  10,  a  detonation  wave  travels 
away  from  the  wall  with  the  speed  6.7  c   =  2.14  km/sec.   It 
raises  the  pressure  to  17.96  atm.  and  the  temperature  to  2980  K. 
The  detonation  wave  is  Immediately  followed  by  a  rarefaction 
wave  through  which  the  pressure  is  lowered  to  4.18  atm.  and  the 
temperature  to  2400°.   The  width  of  the  rarefaction  zone  increases 
proportionally  with  the  time  since  the  beginning  of  the  process. 
Initially,  therefore,  this  drop  is  Instantaneous,  hence,  at  the 
closed  end,  only  the  pressure  4.18  atm.  and  the  temperature  2400 
are  observed.  When   the  detonation  wave  has  reached  the  Interface 


■5J-  The  statement  is  correct  only  inasmuch  as  the  assumption  is 
correct  that  the  detonation  is  a  sharp  discontinuity.   If 
the  finite  width  of  the  reaction  zone  in  a  detonation  were 
taken  into  account  a  continuous  process  would  be  found. 
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chemical  reaction  stops;   a  shock  wave  proceeds  Into  the 
quiet  air  and  a  shock  wave  is  sent  back  into  the  burnt  gas. 
The  shock  in  air  raises  the  pressure  to  18.9  atm.  and  the 
velocity  to  3.16  c   =  1.10  km/sec;  its  ov/n  speed  is  4.04  c   = 
1.41  ton/sec.   The  shock  wave  sent  back  into  the  burnt  gas  is 
extremely  weak.   After  having  traversed  the  rarefaction  wave 
following  the  detonation  front,  this  shock  raises  the  pressure 
from  4.18  to  4.64  atm.  and  the  temperature  from  2400°  to  2430°. 
The  same  happens  at  the  closed  end  as  seen  from  the  diagram  on 
Table  10.   The  rarefaction  wave  that  followed  the  detonation  is 
traversed  by  the  reflected  shock  just  mentioned  and  then  meets 
the  Interface  between  burnt  gas  and  air;   it  is  partly  trans- 
mitted through  the  interface  partly  reflected.   This  reflected 
rarefaction  wave  is  again  reflected  at  the  closed  end,  thus 
decreasing  the  pressure  from  4.64  to  1.93  atm.  and  the  tempera- 
ture from  2430°  to  2070°.   This  reflection  process  takes  place 

between  t  =  .61  t^  and  1.61  t„.   Pressure  and  temperature  at 

o  o  ^ 

the  closed  end  then  remain  constant  up  to  t  =  2.42  t  .when  the 

o* 

next  reflection  arrives. 

In  Mode   DII,    see   Table    11,    the    detonation  wave   travels 
from  the    Interface   towards    the   closed  end.      Its   speed   is   again 
6.17   c      =  2.14  km/sec   and   it  raises   pressure   and   temperature 
to   17.96   atm.    and   2980°K.      A    shock  m.oves    into   the   quiet   air 
raising   the  pressure   to  2.84  atm.    and   accelerating   the  air   to 
a  velocity  of  •815  c    .      Consequently,    the   rarefaction  wave 
which  follows   the   detonation  is   stronger   than  that    occurring    in 
case   DI,    since   it   has    to  accelerate   the  burnt  gas    to  the  velocity 
.815   c      and  to  cause   the   pressure    to  drop   from  17.96   atm.    to 
2.84   atm.      When   the    detonation  arrives    at    the   closed   end,    at 
t   =   .16    t    ,    it    is    reflected   as    a    shock    through  which   the   pressure 
rises    from   17.96   atm.    to  60.90   atm.    and   the    temperature   from 
2980°   to   3675°K.      The   rarefaction  wave    following    the   detonation 
is    traversed  by   the   reflected   shock   and   then    is    reflected  at    the 
closed   end.      During   this    reflection  process  which  takes   place 
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from   t  =   . 16   t      to  t   =   .74   t    ,    the  pressure   at   the   closed 
o  o'  ^ 

end  drops    from  60.90  atm.    (see   Table   3)    to  2.2  atm.    and   the 
temiperature   from  3675°K  to  2210°K.      The   shock  and   the  rare- 
faction wave,    reflected    from   the    closed   end,    will   eventually 
meet   the    interface  between  air  and  burnt  gas.      These  waves   will 
be  partly  reflected,    partly  transmitted.      The   reflected    shock 
is   very  weak;      it   arrives   at   the   closed  end  at   the   time    t  =   1.125 

t     and   is   then  reflected.      The  pressure   at   the   closed  end  has    the 
o 

constant   value   2.2  atm.    during   the    time   interval  from  t   =   .74   t 

^^  o 

to  t  =  1.125  t  .   Due  to  the  reflection  of  the  shock  wave  coming 
in  at  t  =  1.125  t  ,  the  pressure  rises  suddenly  to  the  value 
4.8  atm.,  but  then  drops  quickly  because  of  the  rarefaction  wave 
which  follows  the  shock  wave.   (This  pressure  drop  was  only  esti- 
mated;  it  is  indicated  by  a  dotted  line  in  Tables  1  and  3). 

We  now  proceed  to  describe  in  detail  the  modes  CI  and  CII 
involving  combustion  fronts  with  various  effective  burning  speeds. 

In  mode  CI  combustion  starts  at  the  closed  end  and  pro- 
ceeds into  the  unburnt  gas.   The  combustion  produces  a  compression 
wave  pre -compressing  the  burnt  gas  (see  supplement  [5]  to  the 
Manual).   Since  the  burning  speed  c„  is  asstimed  to  be  constant, 
the  pre -compress Ion  wave  produces  a  constant  compression  and 
hence  is  a  shock  wave.   This  shock  wave  turns  out  to  be  very 
weak.   Nevertheless  it  accelerates  the  gases  to  a  velocity  u 
maich  greater  than  the  burning  speed  c^.   Consequently,  the  speed 
U-  of  the  flame  front  relative  to  the  tube 

U^  =  u  +  c- 

is  also  much  larger  than  the  burning  speed  c_.   The  relationship 
between  the  normal  burning  speed  c„  and  the  flame  speed  U-  will 
be  seen  from  the  Table  in  p.  15.   We  have  in  particular 

U^  =  27.7,   125.   ,  336.  ,  m/sec  . 

Thus,  for  a  burning  velocity  c~  of  3.5  m/sec.  which  may  be  a 
possible  value  for  a  proper  plane  combustion  front,  we  find 
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a  flame-  speed  IT-  =  28  m/aec*   The  flame  speed  TT^  =  336  m/aec 
or  .97  sound  speed  is  rather  high.   Diagrams  for  the  three 
cases  of  mode  CI  are  shown  In  Tables  12,  13,  and  14. 

The  pressure  Increases  through  the  pre-compresslon  shock 
to  the  values  1.10  p  ,  1.53  p  ,  2.88  p  respectively.   The  com- 
bustion process  Involves  only  a  slight  pressure  drop;   It  does 
not  deviate  much  from  a  constant  pressure  burning.   Therefore 
the  pressure  at  the  closed  end  is  still  greater  than  1  atm.  as 
indicated  in  Table  5.   For  the  cases  c-/c   =  .01,  .05,  and  .15 
we  find  p/p  =  1.10,  1.50,  2.52.   Since  the  ccmbustlon  process 
is  slow  the  impulse  keeps  increasing  for  a  long  time;   the 
longer  the  time,  the  smaller  its  rate  of  change  (the  excess- 
pressure).   It  is  not  obvious  in  which  case  the  greatest  final 
value  of  the  Impulse  is  attained,  but  during  the  Intervals  of 
time  considered,  greater  burning  speed  yields  greater  Impulse. 

In  the  case  c  /c   =  .15  the  pressure  of  2.52  atm.  remains  at 
f '  o 

the  closed  end  up  to  the  time  t  =  3.6  t  ;   the  impulse  has 

reached  the  value  I  =  5.51   at  that  time,  see  Table  1. 

o  ' 

In  the  mode  CII  in  which  the  flame  starts  at  the  inter- 
face between  combustible  gas  and  air,  a  shock  wave  is  sent  into 

the  air  accelerating  in  to  velocities  u  =  .03  c  ,  .16  c  ,  .42  c^ 

^^  o'  o'  o 

respectively;      see    the   diagrams    in  tables    12,    13,    and   14.      The 
Interface  between  compressed  air  and  burnt  gas    then  moves   with 
the   same   velocity   in  a   direction  away  from  the   closed  end.      Con- 
sequently,   the   flame   speed  U„  relative   to  the   tube    is    smaller 
than   in  mode   CI,    viz  U„  =   15.2,    76.5,    250  m/sec   as    seen  from   the 
Table   on  p.    15.      The  pre-compresslon  shock  which  proceeds    into 
the  unburnt  gas    is   weaker   than  the   pre-compresslon  shock  for 
mode  CI.      It  raises    the  pressure   only   to  p   =   1.05  p   ,    p   =   1.27   p 
p  =  2.10  p     respectively.      The   pressure   at    the  closed  end,    how- 
ever,   is   greater   than  in  the  mode  CI  because    of   the   reflection 
of  the   pre-compresslon  shock.      The  pressures    there   are   1.10  p   , 
1.60  p^,    4.08   p^  as    shown  in   Table   5.       It    is    interesting   to   see 
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that  the  pressure  at  the  closed  end  for  mode  CII  Is  only 
slightly  greater  than  for  mode  CI  for  c-  up  to  .05  c  .   In 
other  words  the  shock  reflection  at  the  closed  end  essentially 
merely  balances  the  weakening  effect  due  to  the  exhaust  motion 
of  the  burnt  gas  and  compressed  air;   only  for  burning  veloci- 
ties c^  greater  than  ,05  c^  will  the  effect  of  the  reflection 
become  noticeably  dominant.   The  most  remarkable  feature  Is 
that  In  the  case  of  the  high  burning  speed  c -  =  15 c  =  52.05 

/  10 

m/sec   the   Impulse   for  mode  CII   exceeds    that  for   the  detonation 

mode  DI   from  about   t   =  2t       on. 

o 

The  flow  resulting  In  this  case  Is  shown  In  Table  14. 
The  shock  wave  reflected  from  the  closed  end  meets  the  flame 
front  at  t  =  1.01  t^  and  retards  It.   In  addition  to  the  trans- 
mitted shock  wave,  which  Is  very  fast,  a  weak  rarefaction  wave 
is  sent  back  to  the  closed  end  arriving  there  at  t  =  1.23  t  . 
The  pressure  at  the  closed  end  begins  to  drop  slowly  from  that 
time  on.   A  complicated  wave  pattern  ensues  which  Involves  weak 
rarefaction  waves  and  shock  waves  being  reflected  at  the  closed 
end  and  entails  slight  variations  of  the  pressure  there.   When 
the  flame  has  reached  the  closed  end,  which  is  estimated  to  take 
place  at  about  t  =  3  t^,  a  strong  rarefaction  wave  will  start 
at  the  closed  end  and  the  pressure  will  drop  considerably.   It 
is  interesting  to  note  that  the  wave  pattern  here  described  is 
rather  well  confirmed  experimentally  by  Payman's  Schlleren 
photographs,  see  Payman  and  Titman,  [8],  p.  444,  Pig.  4, 

The  discussion  of  the  last  case  of  mode  CII  shows  parti- 
cularly clearly  that  our  way  of  attack  will  in  general  be  success- 
ful  only  for  the  first  stages  of  the  process.   If  the  Information 
for  a  wider  range  of  time  is  desired,  different  methods,  pro- 
bably methods  of  sonic  approximation,  would  have  to  be  employed. 
Since,  however,  the  first  part  of  the  process  yields  the  main 
part  of  the  impulse,  our  Investlgatiftns  appear  to  give  a  signi- 
ficant contribution  to  the  discussion  of  the  effectiveness  of 
certain  jet  operations. 
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